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Abstract 

CN ! We study the initial-boundary value problem resulting from the linearization of the plasma- 

vacuum interface problem in ideal compressible magnetohydrodynamics (MHD). We suppose 
P^ ' that the plasma and the vacuum regions are unbounded domains and the plasma density does 

< 



not go to zero continuously, but jumps. For the basic state upon which we perform linearization 
we find two cases of well-posedness of the "frozen" coefficient problem: the "gas dynamical" 
case and the "purely MHD" case. In the "gas dynamical" case we assume that the jump of 
the normal derivative of the total pressure is always negative. In the "purely MHD" case this 
condition can be violated but the plasma and the vacuum magnetic fields are assumed to be 
Q"^ ' non-zero and non-parallel to each other everywhere on the interface. For this case we prove a 



o 



in 

o 
o 



% 



basic a priori estimate in the anisotropic weighted Sobolev space H^ for the variable coefficient 
problem. 

1 Introduction 

Consider the equations of ideal compressible MHD with the gravitational field ^ G M^: 

dtp + div (pv) = 0, 

dtipv) + dw {pv (^ V - H H) +Vq = pG, 
dtH-V X {vxH) = 0, 
^ dt{pe + ^\H\'^) +dw{{pe + p)v + Hx{vxH)) = p{v,g), 

where p denotes density, ti G R^ plasma velocity, H G W^ magnetic field, p = p{p, S) pressure, 
q = p + ^l-ff P total pressure, S entropy, e = E + ^|wp total energy, and E = E{p, S) internal 
energy. With a state equation of gas, p = p{p, S), and the first principle of thermodynamics, ([1]) 
is a closed system. As the unknown we can fix, for example, the vector U = U{t, x) = {p, v, H, S). 
System ([1]) is supplemented by the divergent constraint 

dwH = (2) 



on the initial data U{0,x) = Uq{x). As is known, taking into account ([2]), we can easily sym- 
metrize system ([T|) by rewriting it in the nonconservative form 



+ divw = 0, 



P 



1 dp 
pc^ dt 

d-fi , ^^ ,_^, ^^ ,. 

(H,V)v + Hdwv 

dt 



dv 
dt 

= 0, 



{H,V)H + Vq = pg, 



(3) 



where c^ = Pp{p, S) is the square of the sound velocity and d/dt = dt + {v, V) (by ( , ) we denote 
the scalar product). Equations ([3]) read as the symmetric quasilinear system 

3 

Ao{U)dtU + Y, MU)djU + Q{U) = 0, (4) 



i=i 



whereQ(C/) = (0,-pe,0), Ao = diag(l/(pc2),p,p,p,l,l,l,l)^ 
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System (JH is symmetric hyperbolic if the the hyperbolicity condition Aq > holds: 

p>0, p,>0. (5) 

Plasma-vacuum interface problems for system ([1]) usually appear in the mathematical mod- 
eling of plasma confinement by magnetic fields. This subject was very popular in the 1950~70's, 
but most of theoretical studies were devoted to finding stability criteria of equilibrium states. 
The typical work in this direction is the classical paper of Bernstein et. al. [2j. At the same 
time, according to our knowledge there are still no well-posedness results for full (non-stationary) 
plasma-vacuum models. Since ([1]) is a system of hyperbolic conservation laws which can pro- 
duce shock waves and other types of strong discontinuities (e.g., current- vortex sheets |14j). it 
is natural to expect obtaining only local-in-time existence theorems. 

The classical plasma-vacuum interface problem models confined plasmas in a closed vessel 
(see, e.g., |lj). In this model the plasma is confined inside a perfectly conducting rigid wall and 
isolated from it by a vacuum region. Let ^~^{t) and ^~{t) are space-time domains occupied by 
the plasma and the vacuum respectively. That is, in the domain i}~^{t) we consider system ([1]) 
(or (21)) governing the motion of an ideal plasma and in the domain 0,~(t) we have the elliptic 
(div-curl) system 

Vxn = 0, div^ = 0, (6) 

describing the vacuum magnetic field "H € M^. Here, as in [21 [3], we consider so-called pre- 
Maxwell dynamics. That is, as usual in nonrelativistic MHD, we neglect the displacement 
current (1/c) dtE, where c is the speed of the light and E is the electric field. 

The boundary of the domain Vt'^{t) is a hypersurface S(t) = {F{t,x) = 0} that is the 
interface between plasma and vacuum. It is to be determined and moves with the velocity of 
plasma particles at the boundary: 

^ = onS(t) (7) 

(for all t € [0,T]). The plasma variable U is connected with the vacuum magnetic field % 
through the relations [21 0] 

[g] = 0, (?^,iV) = 0, {H,N) = ^ onS(t), (8) 

where N = \/F and [q] = q\^ — ^\'H\?j.. These relations together with ([7]) are the boundary 
conditions at the interface S(t). At the perfectly conducting rigid wall T, that is the boundary 
of the vessel fi = 17" (t) U 0,^{t) and the exterior boundary of the vacuum region J7~(t), we have 
the boundary condition 

{n,n)=0 onT, (9) 

where n is a normal vector to T. 

From the mathematical point of view, a natural wish is to find conditions on the initial data 

U{0,x) = Uo{x), xen+{0), F{0,x) = Foix), x E S(0), (10) 



■H{o,x) = noix), xen-{o), (11) 

providing the local-in-time existence and uniqueness of a solution ([/, T-L, F) of problem ([1]), ([6])- 
(jlip in Sobolev spaces. Straightening the unknown interface (see discussion below) and using 
the idea of the partition of unity, this complicated "hyperbolic-elliptic" free boundary problem 
could be splitted, roughly speaking, into two constituents: problem 1^, Q, (fTTj) in the fixed 
domain Cl and problem ([TJ, ©-([H]), p^ . pT]) with an unbounded domain Q~{t). The first 
problem is reduced to the interior Neumann problem for the Laplace equation with a satisfied 
solvability condition by introducing the scalar potential $, where V$ = 7i{t,x). The second 
problem is our main interest in this paper and this problem is a natural generalization to MHD 
the free boundary problem for the compressible Euler equations with the "vacuum" boundary 
condition p\s = (see [6l|T5]). For astrophysical plasmas this problem can be used for modeling 
the motion of a star when magnetic fields are taken into account. 

As in [HI [12], we consider the case of liquid. This means that for problem ([U, ©-([SI), ([TUj) . 
(jlip (with an unbounded domain $7~(t)) the hyperbolicity conditions ([5]) are assumed to be 
satisfied in Q^ up to the boundary S, i.e., the plasma density does not go to zero continuously, 
but jumps. At the same time, in the reality (e.g., for laboratory plasmas [1]) the vacuum region 
is just a region of low enough density. That is, the assumption that the density is small but 
strictly positive at S is quite reasonable. 

Since the interface moves with the velocity of plasma particles at the boundary, at first 
sight the passage to the Lagrangian coordinates to reduce the original problem to that in a 
fixed domain seems most natural. However, as, for example, for contact discontinuities in 
various models of fiuid dynamics (e.g., for current- vortex sheets |14j). this approach seems 
hardly realizable for problem ([T]), ([B])-®, (|10|) . (|lip . Therefore, as in [T3], we will work in the 
Eulerian coordinates and for technical simplicity assume that the space-time domain 0+(t) (the 
plasma region) is also unbounded and the interface S(t) has the form of a graph: xi = ip{t, x'), 
x' = (x2,X3). That is, 

n^{t) = {xi^ip{t,x')} (12) 

and the function Lp{t,x') is to be determined. 

Now we can use a simple straightening of the unknown interface and reduce our problem 
to that in a half-space (see the next section). If, however, the domain il"*'(t) is bounded and 
its initial boundary S(0) is a compact co-dimension- 1 surface in M^, as for shock waves, we can 
follow Majda's arguments [7| (see also [3l sect. 12.4.2]). More precisely, we can make {locally 
in time) a change of variables that sends all boundary locations S(t) to the initial surface S(0). 
We refer the reader to [3 [3] for details of such a change of variables (see also [15] for further 
discussions). 

Thus, we are finally interested in the following free boundary problem. We solve the sym- 
metric hyperbolic system (|3|) (with assumption ([5])) for xi > ip{t,x') and the elliptic system ^ 



[diq] < -e < 0, (13) 



for xi < (p{t, x'). These systems are coupled through the boundary conditions ([7|), ([8]) at the free 
boundary xi = ip{t,x'). Moreover, we have the initial data (jlOp . (jlip (with F = xi — ip{t,x')) 
for t = 0. 

Actually, as for current- vortex sheets [1^, we must regard the last boundary condition in 
([S]) as the restriction on the initial data (|1U|) . More precisely, after straightening of the interface 
and in exactly the same manner as in p^, we can prove that a solution of Q-dH]), (jlOp . (jlip (if 
it exists for all t € [0, T]) satisfies 

divH = in[7+(t) and {H,N) = on S(t) 

for all t G [0, T], if the latter was satisfied at t = 0, i.e., for the initial data (|10p . In particular, the 
fulfillment of div H = implies that systems ([1]) and dH are equivalent on solutions of problem 

o-diD, cnD, dm. 

In the next section we first reduce problem dH^dH])) dlH]), dHI to that in the half-space 
M']_ = {xi > 0, x' gM."^} and then linearize it about a basic state ("unperturbed flow"). For the 
basic state we consider two cases for which we can prove a priori estimates for the linearized 
problem with "frozen" (constant) coefficients. In the first case, for the basic state we require 
the fulfillment of the condition 

dq 

_m_ 

where TV = (—1,0,0) is the outward normal to the boundary of M'^, and [diq] = (diq)\xi=o — 
{v., di'H)\xi=o- Since condition p^ is the counterpart of the natural physical condition dp/dN < 
—e < in gas dynamics (see [6l |T5] and references therein), we call this case "gas dynam- 
ical. " In the second case, we assume that the plasma and the vacuum tangential magnetic 
fields (0,-^2) -f^s) and (0,7^2)^3) are non-zero and non-parallel to each other everywhere on the 
straightened interface: 

|(i/2^3 -^3^2)1x1=0 1 >ei>0. (14) 

We call this case "purely MHD" because for it the physical condition (J13p can be violated, i.e., 
the magnetic field plays a stabilizing role. 

From the mathematical point of view, the principal difference between the above cases is 
that for the "purely MHD" case the symbol associated to the interface is elliptic (see Section 4) 
and for the "gas dynamical" case this symbol can be non-elliptic. We suppose that one can prove 
a local-in-time existence and uniqueness theorem in Sobolev spaces for the original nonlinear 
problem for both of these cases, i.e., "good" initial data for the nonlinear problem reduced to 
that in the half-space Mi^ should satisfy either (J13p or (J14p . In this paper, we manage however 
to prove an a priori estimate for the variable coefficient linearized problem only for the "purely 
MHD" case (see Section 5 for further discussions of open problems). 

The a priori estimate for the variable coefficient linearized problem that we derive for the 
"purely MHD" case can be considered as a first necessary step in proving the local-in-time 



existence for the original nonlinear interface problem by a suitable Nash-Moser-type iteration 
scheme. We plan to use the Nash-Moser method (as in [6l[T5]) because in this a priori estimate 
we have a loss of derivatives. Moreover, the additional difficulty is connected with the fact 
that the interface is a characteristic boundary for the hyperbolic system ^. This implies a 
natural loss of control on derivatives in the normal direction that cannot be compensated in 
MHD (unlike the situation in gas dynamics |1HI15J). Therefore, the natural functional setting is 
provided by the anisotropic weighted Sobolev spaces H^ (see [HI |9l [121 E] and the next section 
for their definition). In this paper we prove our basic a priori estimate in the space Hi. 

The rest of the paper is organized as follows. In Section 2, we obtain the linearized problem 
and formulate main results. In Section 3, for the constant coefficient linearized problem we 
derive a basic a priori L2 estimate for the "gas dynamical" case. In Section 4, for the "purely 
MHD" case we prove an a priori estimate in Hi for the variable coefficient problem. Section 5 
is devoted to the discussion of open problems and contains concluding remarks. 

2 Linearized problem and main results 

2.1 Reduction to a fixed domain 

Let us first rewrite the boundary conditions ([7|), ([8]) for the unbounded domains ()12p : 

dt^ = VN, [q] = 0, nN = 0, Hn = onS(t), (15) 

where vn = {v,N), Hn = {'H,N), Hn = {H,N), N = (1, -52</^, -52<^). 

We straighten the interface S by using the same change of independent variable as in [TJ] 
(see also [15]). That is, the unknowns U and Ti being smooth in Q,^(t) are replaced by the 
vector-functions 

U{t, x) := U{t, $+(t, x),x'), n{t, x) := n{t, $"(*, x),x'), 

which are smooth in the half-space R'^, where 

$^(i, x) := ±xi + *^(t, x), ^=^(t, x) := x(±xi)v?(i, x'). 

and X £ C'^(M) equals to 1 on [—1,1], and ||x'IIloo{k) < 1/2- Here, as in [Ill[l5], we use the 
cut-off function x to avoid assumptions about compact support of the initial data in our (future) 
nonlinear existence theorem. The above change of variable is admissible if 5i$ ^ 0. The latter 
is guaranteed, namely, the inequalities di^~^ > and di^~ < are fulfilled, if we consider 
solutions for which ||vIIloo{[o,t]xR2) < 1. This holds if, without loss of generality, we consider the 
initial data satisfying Hv'oIIlcjoCK^) — ^Z^' ^^d the time T in our existence theorem is sufficiently 
small. 



Dropping for convenience tildes in U and Ti, we reduce (JH), ([6]), (fTSj) . (fTO|) . ([TT]) to the initial 
boundary value problem 

p([/,$+) = 0, V('?^,^-) = in[0,r]xM^, (16) 

]B(C/,'H,c^) = on [0, T] X {xi = 0} X M^ (17) 

iU,n)\t=o = {Uo,no) inRi^, (/^|t=o = '^o inM^, (18) 

where P(C/, ^+) = P{U, '^+)U, 

P{U, *+) = Ao{U)dt + Ii(f/, ^+)ai + A2{U)d2 + ^3(^7)93, 

~ 1 ^ 

^1(^7,^+) = — ^(Ai(C/) - AoiU)dt^+ -Y,MU)dk^^), 

^ k=2 

, , / V xio \ 

5i$^ = ±l + 5i^^, V(?^,^-) = 

y div f) y 
i^ = {nidi^-,nr„nr,), i) = {nn,n2di^-,n3di<^-), 
nn = ni-n2d2^--n3d3^-, -Ur^ = nidi^- + m, 1 = 2,3, 

I dtf -vn \ 



mn,^) 



[q] 
\ nN J 



VN = Vi- V2d2(p - V^d^ip, 



1, 



W\ = Q\xi=o - 2 1^12^1=0' Un = 'Hi- n2d2(p - Hzd-i^p. 
In the MHD equations F([/, '1'+) = we dropped, without loss of generality, the lower-order 
term Q{U) responsible for the gravity. This term is only important for the case of unbounded 
domains for a correct configuration in Sobolev spaces of initial data satisfying conditions like 
(J13p (see ll5j). We just drop this term because it plays no role in our forthcoming linear analysis. 
We also did not include in our problem the equation 

div/i = in[0,r]xM3^, (19) 

where h = {Hn, H2di^^ , Hsdi^'^), Hn = Hi- H2d2^^ - H^ds^^ , and the boundary condition 

Hn = on [0, T] X {xi = 0} x M^ (20) 

because they are just restrictions on the initial data (jlSp (see Section 1). More precisely, referring 
to |14] for the proof, we have the following proposition. 



Proposition 2.1 Let the initial data (|T8|) satisfy (\T9\i and (pO|) . If problem (|T6|) - (|T8|) has a 
solution {U,T-L,f), then this solution satisfies (I19p and (j20p for all t E [0, T]. 

Note that Proposition 12.11 stavs valid if in (J16p we replace system P(C/, ^+) = by system 
([1]) in the straightened variables. This means that these systems are equivalent on solutions of 
our plasma- vacuum interface problem and we may justifiably replace the conservation laws ([T]) 
by their nonconservative form Q. 



2.2 Basic state 

Let 

{U{t,x),n{t,x),^p{t,x') (21) 

be a given sufficiently smooth vector-function with U = (j5, v, H , S) and 

WUWwkUt) + \\H\\wi,(nT) + WdiUWwg^^nT) + Mlwuanr) < ^' (22) 

where i^ > is a constant and 

Qj, ■= (-oo,T] X M3_, af^T := (-oo,T] x {xi = 0} x M^. 

If the basic state (I2ip upon which we shall linearize problem (I16p ~ (ll8p is a solution of this 
problem (its existence should be proved), then it is natural to call it unperturbed flow. The 
trivial example of the unperturbed flow is the constant solution {U,'H,0), where U G M.^ and 
?^ € M^ are constant vectors. 

We assume that the basic state (pTj) satisfies the hyperbolicity condition ^ in 0^, 

p{p, S) > 0, pp{p, S) > 0, (23) 

the first and the third boundary conditions in (I17p on 50^, 

dtif - vn\xi=o = 0, 'Hn\xi=o = 0, (24) 

and system W{7i, ^") = in ilr, 

V X ^ = 0, div ^ = 0, (25) 

where the "hat" values are determined like corresponding values for {U,7i,ip), e.g., 

$^(i, x) = ±xi + $^(t, x), $^(t, x) = x{±xi)ip{t, x'), 

Moreover, without loss of generality we assume that ||'^||Loo(9f2T) ^ ^- This implies 

ai$+ > 1/2, (9i$- < -1/2. 

Note that ([22]) yields 

||T?||H/i(n,) < C{K), 

where W := {U ,diU ,'H,Vt,x^+ ,Vt,x^-), S/t,x = (<9t,V), and C = C{K) > is a constant 
depending on K. 



Remark 2.1 Assumptions (j23P " (j25p are nonlinear constraints on the basic state. We will really 
need them while deriving a priori estimates for the linearized problem. In the forthcoming 
nonlinear analysis we plan to use the Nash-Moser method. As in |141 I15j . the Nash-Moser 
procedure will be not completely standard. Namely, at each nth Nash-Moser iteration step we 
will have to construct an intermediate state (C^„+i/2) 'Wn+i/2) V'n+1/2) satisfying constraints ([25]) - 
(j25p . Without assumption ()25p such an intermediate state can be constructed in exactly the same 
manner as in [HI [15]. Assumption ([25]) does not however cause additional difficulties because 
for given (p it forms together with the last condition in (j24p a boundary value problem reduced 
to the Neumann problem for the Laplace equation in the half-space. We omit corresponding 
arguments and postpone them to the nonlinear analysis. 

Later on, for the linearized problem we will need equations associated to the nonlinear 
constraints (|19p and (j20p . However, to deduce them it is not enough that these constraints are 
satisfied by the basic state (j2ip . As in [14J, we need actually that the equation for H itself 
contained in system P(C/, ^+) = is fulfihed for ([2T]) : 

dtH + —L- I {w, V)H - {h, V)v + ^div u] = 0, (26) 

where 

and w = u — (dt^^ ,0,0). Assume that (I2ip satisfies (I26p . Then, it follows from the proof of 
Proposition 12.11 (see [H]) that constraints (|19p and (|2Up are satisfied for the basic state (|2ip 
if they are true for it at t = 0. That is, without loss of generality we may suppose that (|2ip 
satisfies ([T9]l and (l20D: 

div/i = 0, HnU,=o = 0. (27) 



Thus, for the basic state we require the fulfillment of conditions (l22]) -([2 

2.3 Linearized problem 

The linearized equations for ([16]), ([T7]) read: 

r{U,^+){6U,5^+) := -^F{Ue,^t)\e=o = f in ^t, 

de 

Y'in, ^-)i5n, 6^-) := ^V(?^„ ^-)\s=o = 7- in Qt, 

M'{U,n,0){5U,5n,6^) ■.= ^B{Ue,'He,^e)\e=o = g OnOQT, 

de 

where Us = U + e 6U, Tie = ii + e 6%, c^e = (^ + e 5ip, and 

^f{t, x) := xi±xi)ips{t, x'), ^f{t, x) := ±xi + ^f{t, x). 



6^'^{t,x) := x{±xi)6(p{t,x). 



Here we introduce the source terms / = (/i, . . . , /§), T = {Ti, . . . , J-4), and g = (01,92,93) to 
make the interior equations and the boundary conditions inhomogeneous. 
We compute the exact form of the hnearized equations (below we drop 5): 



"'{U, $+)([/, *+) = P{U, $+)[/ + C{U, $+)[/ - {L{U, $+)^+} 



Y{n, ^-)in, ^-) = Y{n, ^-) + 



diU 



f, 



f VHi X V^- \ 


( " ^ 


(v X 


-n-s 


,v^-) 


\ 


\ ^2 ) 


/ 






i'{u,n,ip){u,n,cp) 



q-{n,n) 

\ T-Ln - 'H2d2(p - Hsd's^ J ,^^ 



where q := p + {H, H), vn '■= vi — V2d2(p — v^dsip, and the matrix C{U, ^^) is determined as 
follows: ^ 

c{u, $+)y = (y, VyAo{u))dtU + (y, VyAi{u, ^+))diU 
+(y, VyA2{u))d2U + (y, Vy^3(f/))93t/, 

, ^ = (yi,---,y8)- 



(y,V,^(f/)):=^,.(^ 



i=l 



Y=U 



The differential operators ¥'{U , ^^) and ¥'(7^, ^ ) are first-order operators in ^+ and ^ 
respectively. As in [TJ], following Alinhac |T], we introduce the "good unknown" 



U :=U 



^H 



5l$^ 



9iC/ 



(28) 



for the hyperbolic system of linearized MHD equations. Similarly, we also introduce the "good 
unknown" 



'H:=n 



^ 



di^- 



■diH 



(29) 



for the elliptic system for the perturbation of the vacuum magnetic field. Taking into account 
assumptions (p^ and ([25]) and omitting detailed calculations, we rewrite our linearized equations 
in terms of the new unknowns (12811 and (1291) : 



P([/,^+)f7 + C(C/,^+)[/ 



^H 



Yin,^-) = T. 



ai{L(c/, $+)}=/, 



(30) 
(31) 



10 



W{u, n, ^){u, n, if) := n'iu, n, <^)(c/, n, ^) 

' dt^ + 7)292'/' + v^d?.'^ -VN - 'fdiilN 
q - i^n) + [diq]^ 

where vn = vi — 1)292'^ — vsdsip, Hn = vi — 'H2d2'f> — 'H?.,dz0, and 



5, (32) 



[diq] = (d: 



iq)\xi=o 



i'H,din)\,,=o. 



We used the last equation in (j25p taken at xi = while writing down the last boundary condition 
in dMD- 

As in [U [T^ 115] . we drop the zeroth-order term in ^+ in (|30|) and consider the effective 
linear operators 

P;({7, ^+)U := P{U, ^+)U + C{U, $+)[/ = /. 

In the future nonlinear analysis the dropped term in (|30p should be considered as an error term 
at each Nash-Moser iteration step. 

Regarding system (|31|) . without loss of generality we may actually drop the source term T. 
At first sight, we have to keep it because the nonlinear system Y{7i, ^^) = will produce errors 
in the Nash-Moser iteration scheme. That is, in the future nonlinear analysis we will have to 
go outside the class of divergence-free irrotational fields. At the same time, it follows from the 
detailed analysis of an exact form of the accumulated errors for the elliptic system Y{l-L, ^^) = 
and the boundary condition 'Hn\xi=o = (corresponding arguments are omitted and postponed 
to the nonlinear analysis) that the source terms J- and (73 have the following special form: 



T 



V X «B 
divb 



93 — ^l|x-i=05 



(33) 



where 



<3 = ihdl<^-,br„br,), b = ibn,b2dl<^-,hdl<l>- 



bodo'i' - 63 ^3^ 



bidi^ +bi, 



2,3, 



and b{t,x) = (^1,621^3) is a vector-function. Passing to the new unknown 

n' = 'H-b 
and omitting then the primes, in view of (j33p . we get the homogeneous system 

V(-H,$-) = 
and the last boundary condition in (|32p becomes homogeneous (53 = 0). 



11 



We now write down the final form of our linearized problem for {U,'H,<p): 

3 

AodtU + Y^ Ajdjil + CU = f in Qt, (34) 



(35) 



dt(f = VN - V2d2(p - V'id'i'^ + fdiVN + gi, 

q = {'H,'H) - [diq]ip + g2 on dQr, 

V X i^ = 0, div f) = in Ot, (36) 

-Hn = d2 {n2v) + ds {-Hs^) on d^r, (37) 

{U,'H,ip) = fort<0, (38) 



where 



A^=:A^{U), Q = 0,2,3, Ai =: Ai{U,^+), C:=C{U,^+), 

■Hn = Til - n2d2^- - ■H3d3^- , 'Hr.^=nidS-+'Hi, i = 2,3. 

We assume that / and g = (51,52) vanish in the past and consider the case of zero initial data, 
which is the usual assumption. We postpone the case of nonzero initial data to the nonlinear 
analysis (construction of a so-called approximate solution) . 

2.4 Basic a priori estimates 

We first write down our basic a priori estimates for the case of constant ( "frozen" ) coefficients 
of problem (|34|) - (j38p . Before formulating this result we give the definition of the anisotropic 
weighted Sobolev spaces i7™. Following [TBI 13 HI E]) the functional space H^ is defined as 
follows: 

iJf (r3_) := L G L2{M.l) I d^d^u G L2(M^) if \a\+2k<m\, 

where m G N, 9" = {adi)°'^d2^d^^ , and o"(xi) G C°°(M-|-) is a monotone increasing function 
such that cr(xi) = j;i in a neighborhood of the origin and cr(xi) = 1 for xi large enough. The 
space -ff™(M5|_) is normed by 



m,* - 2^ \\(^*'^i'"'\\l2{RI)- 
\a\+2k<m 



We also define the space 

m 
k=0 

equipped with the norm 

„T rn 

bAm,*,T= WHtm'i^^dt, where |||u(t)|||^., = J]] ||a^u(t)|| 

•^-°° i=o 
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2 

m—j,*' 



Within this paper we use the space H^{Qt) mainly for m = 1. Clearly, the norm for HI{Qt) 
reads 

[u]l*,T = [ {u^ + {dtuf + {adiuf + {d2uf + {dzuf) dtdx. 

We are now in a position to state our main results. 

Theorem 2.1 Let the basic state (|21|) satisfies assumptions (|23p - (|27|) . Let the coefficients 
of problem (J34p - (j38p are "frozen", i.e., the coefficients of the interior equations ()34p . (|36p 
and the coefficients of the boundary conditions (|35p . (j37p /lawe 6een calculated at given points 
(t*,j;*,x'^) G ilr o-nd (t*,x'^) € 9^7^ respectively, in particular, the coefficient [diq] is a constant. 
Then, for the "gas dynamical" case (J13p , 

[dig] > 0, (39) 

sufficiently smooth solutions {U,T-l,(p) of (jM|) - (j38P o6e|/ the estimate 

WllLzinT) + II^IIl2(!^t) + Wv^WiiidnT) ^ C"! 11/11^2(^7) + WaWmidnr)} > (^^0) 

where C = C{T) > is a constant independent of the data {f,g)- 

For variable coefficients, we have not yet managed to derive an a priori estimate for the "gas 
dynamical" case, but for the "purely MHD" case we prove the following theorem. 



Theorem 2.2 Let the basic state (j2ip satisfies assumptions (J22p - (j27p . Then, for the "purely 
MHD" case <^ . 

|(#2?^3-i^3^2)Ui=o|>ei>0, (41) 

sufficiently smooth solutions {U,T-L,(p) of problem (J34p - (j38p obey the estimate 

[U]i,*,T+ WHWminr) + \W\\m{dnT) ^ C'{[/]2,*,t + 1151^2(9!^^,)}, (42) 

where C = C{K,T) > is a constant independent of the data {f,g). 

3 "Gas dynamical" case 

3.5 Properties of problem (j3il) - ([38l) 

Before "freezing" coefficients we discuss some useful properties of the variable coefficient problem 
P4p -()38 p . First of all, as for current-vortex sheets |14] . we can prove the following proposition 
(see Appendix A in [14J for the proof). 

Proposition 3.1 Let the basic state (j2ip satisfies assumptions ()22p -()27 p . Then solutions of 
problem ([MD^dMD satisfy 

div/i = r in^T, (43) 
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H2d2(p + H^d^ip - Hn - f diHjsf = g^ on dVir- 



(44) 



Here 

h = {Hn, H2di^+, H;dil+), Hn = Hi- i?292$+ - H^di^+ 

(Hi\l\x-^=Q = Hn\xi=o), the functions r = r(t,x) and 53 = g^{t,x'), which vanish in the past, 
are determined by the source terms and the basic state as solutions to the linear inhoniogeneous 
equations 

(45) 



1 



dta -\ ^^ {{w,\/a) + adivu} = T inVt^^ 



5*93 + V2d2gz + hdzgz + (92^2 + ^3^3) gi. = G on OQt, 
where a = r / di^^ , T = {div f h) / di^^ , 

fu = ifn, /6, /r), /n = /5 - hd2^^ " fjO^^^ , Q = {d2{H2gi) + 93(^35l) " /«} |,,=o 



(46) 



It follows from the first condition in ()24p that the interior equation ()45p does not need a 
boundary condition because wi|xi=o = 0. Therefore, from (j45|) we get 

lk||L,(Q,) < CWfWH^^n,) < C[f]2,.,T. (47) 

Here an later on C is a constant that can change from line to line, and sometimes we show the 
dependence of C from another constants. In particular, in (|47p the constant C depends on K 
and T. Using (|46p and the trace theorem [9] for the spaces H^, we easily estimate: 



15311^1(903,) < C" {1151 11^2 (an-r) + \\f\xi=o\\m{dQT)} - ^ihWH^idUT) + [/]2,*,t}- (48) 



In view of the first condition in (124p and the second condition in (I27p , the boundary matrix 
Ai is singular at xi = 0, i.e., the plane xi = is a characteristic boundary for the hyperbolic 
system (j34p (exactly as for current-vortex sheets |14]). Following |14] . we introduce the new 
unknown 

V = {q,Vn,V2,V3, Hn, H2, H3, S) 

for separating "characteristic" and "noncharacteristic" unknowns. We have U = JV, with 



J 
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where Hj-^ (= 2, 3) are determined in the same way as 71^ above. Then, system (plj) is equiva- 
lently rewritten as 

3 

AodtV + J2^kdkV + A4V = F, (49) 

fc=i 

where Aa = J''^ AaJ {a = M), F = j"""/, and 



^4 = JT J c + AodtJ + ^ AAJ I 



The boundary matrix Ai in system ()49p has the form 

^1 = ^ + ^(0), ^=^^^12, ^(o)Ui=o = 0, (50) 

where £ij is the symmetric matrix which (ij)th and {ji)th. elements equal to 1 and others are 
zero. The explicit form of A(q\ is of no interest, and it is only important that ^(o)|xi=o = 0- 
Therefore, the boundary matrix Ai on the boundary xi = is of constant rank 2. That is, (j49p 
is a symmetric hyperbolic system with characteristic boundary of constant multiplicity (in the 
sense of Ranch [10] ) . It is also noteworthy that because of (j43]) not only q and Vn but also Hn 
is a "noncharacteristic" unknown. For the "noncharacteristic" part of the vector V, 

Vn = iq,Vn,Hn), (51) 

we expect to have a better control on the normal (xi-) derivatives. 

We now discuss the elliptic part of our problem. In [2], the vacuum vector potential was used 
for the div-curl system. Unlike [2], here we introduce the scalar potential A for system l\36\i : 

S^ = VA. (52) 

Then, we get 

AA = 0, (53) 

with ^ 

A = df + di + dl 5i = -^ai, dk = dk-^^di, A; = 2,3. 

Passing to the "original" curvilinear coordinates xi = ^~{t,x), x' = x', we could rewrite (f53|) 
and the boundary condition (j37p as the Neumann problem for the Laplace equation for A (with 
a satisfied solvability condition). However, we do not need do so. Moreover, we do not rewrite 
(j37p in terms of the potential A. In fact we will only use relation (j52p and it will be even more 
convenient for us to work with the equation 

div ^ = (54) 

instead of equation ([53]) . 
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3.6 Proof of Theorem \2A] 

We first do not "freeze" coefficients and obtain an inequality for variable coefficients which will 
imply the a priori estimate (j40p for "frozen" coefficients. By standard argument we obtain for 
the hyperbolic system (|49p the energy inequality 

ft 



(55) 



I(t)-2^^ '?%U=odx'ds<c(^||/||i^(^^) + ^ I{s)dsy 
where I{t) = Jj^s (^o^i V) dx. It follows from the boundary conditions ([35]) that 

-2qvN\^^=o = 2 {[diql^p - 92) vn\xi=o - 2{'H,'H)vn\xi=o, (56) 

where 

2{[diq\ip- g2)vN\xi=o 

= 2{[diq]ip - 52)(«9t(/3 + i)292V3 + vs^ac^ -ipdiVN - fi'i)U'i=o 
= dt{[diq]ip'^ - 2g2(p} + 82 {v2\xi=Q [diq]'^'^ - 2^2 1 a; 1=0 529'} 

+ d'i{v'i\^^=Q[diq]ip^ - 2v^\^^=Q g2ip] +2gig2 
- {dt[diq] +d2{v2[diq\) + d^{v^[diq]) - 2[diq]div n]\x^=q V^ 

+ 2 {dtg2 + d2{v2g2) + d3{v3g2) + g2diVN - [diq]gi} Ui=o ^- (57) 

Assume that [diq] > e > (this is the version of condition ([39]) for variable coefficients). 
Using the Young inequality, from (j55|) - ([57|) we obtain 



Ht) + l: [ [diq\ipUx'-2 [ {'H,n)m\x^=odx'ds 

< C{K,e){\\f\\l^^^) + ||<7llii(an.) + / {ns) + ||V'(^)||L(M2)) ds}. (58) 



Regarding the boundary term —2{'H,'H)vn\xi=o in ([58]) . in view of ([35]) . ([37]) . ([52]) . and the 
second condition in ([24p . we have 



- 2{'H,'H)vn\xi=o = -2{'H2d2 A + nads A) {dtif + V2d2^ + hdsy:} - (f diVN - gi)\xi=o- (59) 

Note that in the framework of the L2 theory we are not able to treat the term 2gi (Ti, %) \xx=o 
contained in ([59p directly. On the other hand, since in (|58p we anyway lose one derivative from 
g, we can use the classical argument suggesting to reduce our problem to one with homogeneous 
boundary conditions by subtracting from the solution a more regular function. Namely, there 
exists U = {p,v,H,S) G H^{Qt) (or more precisely, {q,Vn) G H^{Qt) and {v2,vs,H, S) G 
HI{Qt)) vanishing in the past such that 

-VN = 91, q = 92 ond^T, 
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where {;„, vn, and q are determined like corresponding values for U. If [/ = [/^ + U, then 
{U^n, ip) satisfies ([MD-dMD with 5 = and / = f\ where f^ = f- ¥'^0, $)[/. That is, it is 
enough to prove estimate (|40p with g = 0. Without loss of generality, we will just assume that 

in dMD 91 = 0. 

Let us consider the term —2dtf{7i, 'H)\xi=o contained in ([M|) . Integrating by parts and using 
(1371) and (1521). we obtain 



-2/ {n2d2A + n3dsA)u,=odt^dx'ds 

JdUt 

= -2 [ (VA, {dti^n') - ^dtU') 

Jant ^ 

= 2 [ AdtnNU=odx'ds + 2 [ {S^',dtn%,=oipdx'ds, (60) 

JdUt JdUt 

where V' = (82,83), %' = ('^2,^3), and S^' = {'Hr2,T^Tz)- It is noteworthy that for "frozen" 
coefficients the last integral in ([60j) disappears. Regarding the penultimate integral in ([60]) . 
taking into account (j36p and (j52p . we have 



2/ AdtnN\x,=odx'ds = J{t)-2J ■HNdtA\^,=odx' ds 

JdUt JdUt 

= J{t) + 2 I di {nndtA) dx ds 
= J{t) + 2 f ^-Hndtinidi^') -dtAdiv^,{n'di^-)^dxds 
= J{t) + 2 f {-Hndti-Hidi^^) 
+ 9i$- {n', dtCHiV^- + n')) } dx ds = J{t) - K{t) + L{t), (61) 
where H' = (7^2, ^3), div^./6' := 9262 + ^363 (with h' = (62, h) G K^), 

J{t) = 2 [ AnN\x^=odx', K{t) = - [ di$~\'H\'^dx, 

L{t)= I \didt^-{nl-\n'\^) + {{di^-Vdt^- -didt^-V^-),-n')^dxds. 

Recall that di%- < -1/2, i.e., K{t) > i \\'H{t)\\l,^3 y For "frozen" coefficients L{t) = 0, 
but even for the variable coefficients case we easily estimate —L{t) from above by C||7t(|[^ ,^ - 
Multiplying equation (j54p by the potential A, integrating the result over the domain M^, and 
using then integration by parts and ([52|) . we get 

J{t) = 2K{t). (62) 

Thus, it fohows from dnUD-dS^D that 

- 2 / dt^{n,n)\x,=odx' ds > K{t) - C||7^|li^(^^) > I II^WIIL(R3) - Cml^^^y (63) 
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Consider now the term — 2({)', VV)(^)^)U'i=o (with v' = {v2,'iJ3)) contained in (j59]) : 



2 / {n2d2A + n3dsA){v',V^)U,=odx'ds 

JdUt 



= -2 [ \{v',VA){n',v^) 

+ (hn^ - V3'H2){d3Ad2ip - d2Ad3^)} dx' ds = M{t) + N{t), (64) 



where 



M(i) = -2 / ■HN{v\mx,=Gdx' ds, 
JdUt 



N{t) 



an, 



{ (93(^)2^3 - ^3^2) + *25i(?^jv)) -Hr, 

+ (92(1)2^3 - i)3^2) + ^3^1 (^7V)) Tir,} 



xi=0 



ipdx' ds. 



For "frozen" coefficients N(t) = 0. To avoid unnecessary technical details we consider here the 
integral M(t) for the particular case ip = and leave the general case to the reader. For ip = 

M{t) = -2 [ {V2nin2 + *3^1^3) Ui=0 dx' ds. 

Jaut 

and the div-curl system (j36|) takes the form V x 'H = 0, dwH = 0, with 7i = {—'Hi,'H2,'H3)- 
From this system we easily deduce that 



Jai 



'^i?T(2Ui=odx'ds = and 



'Hi'H3\xi=o dx' ds = 



I ant JdUt 

For the "frozen" coefficients case, V2 and V3 are constants and we therefore conclude that M{t) = 
0. For variable coefficients, omitting simple calculations, from the div-curl system we derive the 
estimate 

M{t) > -Cml^^^y (65) 

At last, consider the term 2{p{diV]\f){'H, 1-1)1x1=0 contained in ([59]) . For variable coefficients, 
in view of assumptions (|24p and (j25p . the corresponding boundary integral together with N{t) 
and the last integral in (I60p can be written in the following compact form: 



N{t) + 2 / {{n',VA)diVN + {^i',dtn')} 



a;i=0 



ipdx' ds 



N{t) = 2f (S)',{dt'H-V xivxii)}') ^dx'ds, (66) 



where {6}' = (62,63) for any three-dimensional vector h = (61,62163)- Since we cannot control 
the trace 'H.\x^=o we do not know yet how to estimate the integral M(t) (see Section 5 for further 
discussions). But for "frozen" coefficients, taking into account ()62p . we have 

M(t) = 2{diVN) I {n', VA) 1,^=0 ^ dx' ds = -2(^1%) / K{s) ds. (67) 

Jant Jo 
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It follows from (EHIl-dMD that 



m + \ P(i)llL(M^) + \ fjdiq] ^2 dx'+M{t) 

+ f {i{s) + ii?:f(.)iii,(i,3^) + ii¥'(^)iiL(R^)) ds}- 



'0 

For the "frozen" coefficients case, by virtue of (j67p . we obtain inequality (|68p where we formally 
set M{t) = 0. Taking into account assumption (j39p and applying Gronwall's lemma, from 
this inequality we finally deduce the basic a priori estimate (|4Up . This completes the proof of 
Theorem [O 

4 "Purely MHD" case 

For the "purely MHD" case thanks to assumption ()4ip we can resolve ()35p . (|37p. and (j44p for 

^t,x'^ = (5*9?, VV): 

Vi,^/(/9 = aiii'jv|:ri=o + «2^Ar|xi=o + a3^'Ar|xi=o + a4¥' + ao53, (69) 



where the vector-functions a^ = a«(W^|xi=o) = (^a)^a'^a) can be easily written down, in 

3 — "3 



particular, ai = a\ = 0, 



a-, = ^ ^ ^^^^ , Oo = ^^^^ ^^^^ , etc. 

{H2'H3 - H^'H2)\x,=0 {H2n3 - i?3^2)Ui=0 

Using the terminology of paradifferential calculus, we can say that for the "purely MHD" case 
the symbol associated to the interface is elliptic (see, e.g., [3j). This fact plays the crucial role 
in the proof of estimate (j42p . 

Using the special structure of the boundary matrix Ai (see ([50]) ). from system ([I9|) . equation 
(jl3|) . and estimate ([I7|) we easily deduce the inequality 

WdiVnimU^l) < C{K){[f]l,^T + |||f7(t)|||t} • (70) 

Moreover, by resolving the div-curl system ()36p for the normal normal (xi-) derivatives, one gets 

wdinmi^^s^^ < c{K){\\d2m\\l^^i) + ii53^(t)iiL(M3^)} • (71) 

Then, thanks to the trace theorem ()70p and (j7ip imply 

llKU=o(t)|lL(M2)<c{[/]|,,T' + ir(t)lllv}, (72) 

||H„U,=o(t)|lL(M.) <c{\\d2m\\U^l^ + l|53?i(0llL(R^)} • (73) 



19 



We first get an estimate for weighted derivatives. To estimate such terms we do not need 
boundary conditions because the weight cr\x^=o = 0. By applying to system ([M|) the operator 
adi and using standard arguments of the energy method (see, e.g, [13j for more details), we 
obtain the inequality 

lk^if>WllL(Mp <G{[/]t,T + [f>]?,M}- (74) 

We can easily get the inequalities 



and 

where (|75|) follows from the trivial relations 

d ..^- ,2 _o f fTT ^T■r^A^ d ,, .,,,,2 



di WUmU^i) =2 1 iU,dtU)dx, - ||W(t)||t(j,3^) = 2 1 in,m)dx, 

and (I76p is the result of the multiplication of the first boundary condition in (I35p by 2^3, the 
integration over the domain M^, and the usage of ([69]) . estimates ([72]) and (jiHj) . In ([76|) we use 
the norm 

llk(0lllHi(M2) := / (n' + {dtnf + (92^)^ + {d,,uf) dx' . 

We now proceed to estimating the tangential derivatives 9^ and dt {k = 2, 3) of the solution. 
This is the most important step because we shall use the boundary conditions. Differentiating 
system pUj) with respect to Xi (with £ = 2 or i = 3) and using again standard arguments, we 
get the energy inequality 

f {AodeV,deV)dx-2 f deqdiVN\x,=odx' ds < C \[f]l,T + [U]Lt} , (77) 



By using the boundary conditions (|35p we obtain 

-2/ diqdiVN\x,=odx'ds = Ji{t) + J2it), (78) 

JdUt 



I an. 
where 



dx'ds, 

a;i=0 



Ji(t) = 2 / \{[diq]deip-deg2)deVN + [dediq]ipdeVN - {de'H,'H)deVN\ 

J2{t) = -2 / {ii,dm)diVN\x,=odx' ds. 
J ant 

Integrating by parts and applying (j69|) . we get 

-^1 {t)= I \ ciVNdiip + C2VNgz + csVNdegs + 2vNdjg2 
Jaut ^ 

+ c^HNdgVN + y^ Ci+iHNdiVN + csvn^ (■ dx' ds, 

^-^i=l J 1x1=0 
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where 

Cj = Cj{W\^^=o, diW\^^=Q) (j =TJ), eg = cs(W\^,=o, diW\^^=Q, a|t?|^.^=o) 

are functions (coefficients) dependent on the basic state (pT|) . To treat the terms c^HNd^VN and 
Ci+iHNd(,VN contained in the boundary integral Ji(t) we use the same arguments as in p^lTl] . 
That is, we pass to the volume integral and integrate by parts. In particular, we have: 



JdUt Jut 

= / \ c^deHNdiVN + {deC5)HNdiVN 

- c^diHNdeVN - idiC5)HNdiVNjdxds, (79) 

where C5\xi=o = £5. Taking into account ([70]) - ([72|) and ([l8j) . we can now estimate the boundary 
integral Ji(t) as follows: 

Regarding the integral J2{t), to treat it we can repeat arguments in ([59]) - ([66|) , where, roughly 
speaking, instead of A, 93, and Tij^ we have d^A, dup, and diV-N and there appear additional 
lower-order terms. In particular, the integral 

/ {coeS\di^')l^^^d,ipdx'ds (81) 

JdUt 

is the counterpart of the integral M{t) in ([66]) . where coeff' = (coeff2, coeffa) and coeffj are 
coefficients dependent on the basic state (|21|) . But now we can treat such an integral for variable 
coefficients by using the ellipticity of the interface symbol. That is, we express dgip appearing 
in ([HT]) through Hj^^^^^q, 'Hj^\xi=q^ ^^^d §3 (see ([69]) ). After that we use the same arguments as 
in (j79p . Omitting detailed calculations (in particular, we also exploit (j73p ). we finally estimate 
the boundary integral J2(t): 

+ [U]Ut + \\njj.(^n,) + Mm(dn,)}- (82) 
It follows from ([77]), ([TSD, ([80]), and (l82|) that 

+ [U]Ut + m mm + Mmidn.)}- (83) 
Combining ([83]) with i^, dHD-dH]), we get 

< C(^){[/]t,T + Mh^anr) + P]U + miming + Mmian,)}- (84) 
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We still miss the L2 norms of dtU and dfTi in the left-hand side of (j84p which we need for 
closing estimate (|84p. To estimate the time derivatives of U and "H we use the same arguments 
as in ([77 |) " ([83]) . We just replace di by dt (or dg), and the only principal difference is that ex- 
pressions in the form dt{- ■ ■) do not disappear after the integration over the domain Qt (whereas 
/q di{- • • ) dxds = 0). For example, the right-hand side of ([7^ with df replaced by dg contains 
the additional integral 



/ ds{c5HNdiV]\f)<ixds = — / c^HNdiVN dx. 



'fit 
Using the Young inequality, (|70p , and (j75p . we estimate this integral as follows: 

- / c,HNd,VNdx<CiK)\[f]l,^T + e\\\Um\l. + -[U]l,,t], 

where e is a small positive constant. Dealing analogously with remaining terms appearing after 
the integration by parts with respect to t (see also [13] for similar calculations), we finally obtain 
the inequality 

mml^^i) + imml^^i) < c{K){[f]i^^T+\\9\\i^9nr) 

I r7'7']2 I \\nj ||2 i 11 ,^||2 

+ ['-'ii,*,t + \\'~''\\m(Qt) + IWWmidQt) 

+ e||[/(t)||?_, + e||?i(t)||^,(jj^)}. (85) 

Combining (j85p with (j84p and choosing the constant e to be small enough, we get 



\t + iiimm^H^ui) + \MtWm 

< ^(i^){[/]t,T + Mnnan^) + [U]l*,t + PII?^i(no + Mmiant)}^ (86) 

where 

\M*)\\\m(Rl) ■= \H't)\\m(Ri) + ll^t^(*)llL(Mp- 

Applying Gronwall's lemma, from the last inequality we derive the desired a priori estimate 
This completes the proof of Theorem 12.21 



5 Concluding remarks/Open problems 

We have obtained first results towards the proof of the local-in-time existence of smooth solutions 
of the plasma- vacuum interface problem in ideal compressible MHD under the basic assumption 
that the plasma density does not go to zero continuously, but jumps. We have found two cases 
of well-posedness of the constant coefficient linearized problem: the "gas dynamical" case and 
the "purely MHD" case. From the mathematical point of view, for the "gas dynamical" case 
the interface symbol can be non-elliptic and for the "purely MHD" case the interface symbol is 
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always elliptic. For the latter case the MHD counterpart (|13|) of the natural physical condition 
in gas dynamics [H [15] can be violated, i.e., the magnetic field plays a stabilizing role for well- 
posedness. 

For the "purely MHD" case we have managed to derive a basic a priori estimate for the 
variable coefficient linearized problem. We prove this estimate in the anisotropic weighted 
Sobolev space Hi because the interface is a characteristic boundary for the hyperbolic system 
of MHD equations, and in MHD a natural loss of derivatives in the normal direction in a priori 
estimates cannot be compensated as in gas dynamics [15j. Assuming that the original nonlinear 
problem has smooth enough solutions, we can easily prove the uniqueness of a solution of this 
problem by a standard argument and using the basic a priori estimate in Hi. 

In the basic a priori estimate (j42p we have a loss of derivatives from the source terms to 
the solution. Clearly, we will have a loss of derivatives also in a corresponding tame estimate 
whose derivation is postponed to the future. Therefore, we expect to prove the existence of 
solutions of the nonlinear problem by a suitable Nash-Moser-type iteration scheme. We do not 
see any obstacles in this direction. We think that the derivation of the tame estimate and the 
realization of the Nash-Moser procedure is just a technical matter and can be done as in |14j for 
current-vortex sheets. 

At the same time, there is still an open problem in getting the local-in-time existence result 
for the "purely MHD" case. The point is that we have not yet proved the existence of solutions 
of the linearized problem. We can naturally formulate a dual problem for it, but we still do not 
know how to get an a priori estimate for the dual problem. This is a very surprising fact because 
usually if we can obtain an a priori estimate for the original linearized problem, then in exactly 
the same manner we can derive it for the dual problem. After that the existence of solution of the 
linearized problem can be proved by the classical argument of Lax and Phillips [5j. Of course, 
our "hyperbolic-elliptic" problem is very nonstandard and this causes the mentioned difficulty. 
We expect to prove the existence of solutions of the linearized problem either by iterations or 
by considering a regularized problem. This work is postponed to the future. 

Regarding the "gas dynamical" case, it is still unclear how to carry the basic a priori estimate 
obtained for the constant coefficients case over variable coefficients. The difficulty is connected 
with the appearance of additional lower-order terms and the fact that we cannot control the trace 
of the perturbation of the vacuum magnetic field in the higher norm. We think that it is very 
unlikely that the plasma-vacuum interface problem is not well-posed for the "gas dynamical" 
case. But, we cannot completely exclude this possibility. This question is the most important 
open problem both from the mathematical and the physical points of view. On the other hand, 
for the m,odel free boundary problem when the vacuum magnetic field H = we can prove a 
local- in-time existence theorem in the anisotropic weighted Sobolev spaces H'^, provided the 
initial data satisfy the condition dq/dM < — e < (cf. (fT3]) ). Roughly speaking, we can prove 
such a theorem by the combination of arguments applied in [15] to the compressible fluid- vacuum 
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problem and in |14j to ideal compressible current-vortex sheets. 
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